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Spin-polarization is known to lead to important magnetic and optical effects in open-shell atoms 
and elemental solids, but has rarely been implicated in controlling structural selectivity in com- 
pounds and alloys. Here we show that spin-polarized electronic structure calculations are crucial for 
predicting the correct T = crystal structures for P&zX and Pt$X compounds. Spin-polarization 
leads to (i) stabilization of the LI2 structure over the DO22 structure in PtaCr, PdaCr, and PdsMn, 
(ii) to the stabilization of the DO22 structure over the LI2 structure in PdaCo and to (iii) order- 
ing (rather than phase-separation) in PtsCo and PdaCr. The results are analyzed in terms of 
first-principles local spin density calculations. 

PACS numbers: 61.66.Dk, 71.20Cf, and 75.50.Cc 
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Crystal structure compilations |l]j2|] reveal that the 
most commonly occurring structures among intermctal- 
lic binary compounds with a 3:1 stoichiometry (A3B) 
are the cubic LI2 and the tetragonal DO22 (Fig. |l|). The 
crystallographic difference between the LI2 and DO22 
structures is rather subtle: the two structures have iden- 
tical first neighbor coordination (each A has 8A + AB 
neighbors and each B has 12A neighbors) while a differ- 
ence exists in the second shell (see Fig. |l|). The manner 
in which particular A3B compounds select the LI2 or 
the DO22 configuration appears to be rather interesting. 
For example, || the 4c? trialuminides AI3M show the 
sequence LI2 — * LI2 — » DO22 as M varies across the 
4d row Y— s-Zr— >Nb, while the 3c? palladium alloys Pd^X 
show LI2 — ► LI2 — > DO22 — * LI2 — » LI2 — ► LI2 as one 
proceeds in the 3d row X = Sc— »Ti— >V— »Cr— »Mn -^Fc 
|| (for X — Co and Ni, the systems phase-separate). The 
origin of such regularities was the subject of numerous 
investigations including the c?-electron "generalized per- 
turbation method" (GPM), [Q and first-principles calcu- 
lations. (^-0] However, these calculations failed to repro- 
duce the observed structural trends. These calculations 
were non-magnetic (NM, i.e., without spin-polarization). 
This appeared to be a reasonable assumption, since one 
expects that an alloy rich in a non-magnetic compo- 
nent (e.g, PdsCr) or one without any magnetic com- 
ponents (e.g, Pd3 V) will not have any significant mag- 
netic effects. We demonstrate here that spin-polarization 
have a crucial influence on the structural stability of 
Pd$X and Pt$X compounds: it stabilizes the observed 
LI2 structure over the DO22 structure in PtaCr, PdsCr, 
and PdaMn, the DO22 structure over the LI2 structure 
in PdaCo, and is responsible for compound formation 
(rather than phase-separation) in PtaCo and PdaCr. 

The key insight to stability in compounds and alloys 
has traditionally been the association of stability with 
low density of states (DOS) at Fermi energy Ep. M 
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FIG. 1. The crystal structures of the (a) LI2 and (b) DO22 
structures. The insert shows the atomic coordination about 
A and B sites in the first (1st) and second (2nd) atomic shells. 
Ai and A2 indicate two distinct A sites in the DO22 structure. 



Nicholson et al. |p[ noted that for transition metal alu- 
minides the more stable of the two structures (LI2 or 
DO22) corresponds to the one with smaller DOS at the 
Fermi-level, N(Ep). Local density approximation (LDA) 
H band structure and total energy calculations ||]7|] have 
later substantiated this relation. To examine such a re- 
lation for inter-transition metal A3B compounds rather 
than aluminides we have calculated N{E) (Fig. ||) and to- 
tal energy difference SE = E(L1 2 ) - E(D0 2 2) [Fig. |(a) 
and Table § for Pd 3 A with 3c! atom X = Sc through Cu 
using the NM linearized augmented plane wave (LAPW) 
method. |ll| We see that the structure with the lower 
calculated NM total energy (Table |) indeed has a lower 
calculated NM N(Ep) (Fig. 2), thus substantiating ear- 
lier trends for aluminides. P - ?} For example, for Pd$X 
with X = V, Cr, and Mn, the Fermi energy, Ep, falls 
near a DOS maximum for LI2 but near a DOS mini- 
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mum for DO22] correspondingly E{D022) is lower than 
E{L\2). Unfortunately, while the magnitude of N(Ep) 
is indicative of the stability of the calculated structure, 
these non-magnetic calculations incorrectly predict the 
observed stable crystal structure in several cases: while 
Pd 3 V is correctly predicted to be more stable in the DO22 
structure, the observed stable structure for Pd 3 Cr, 
Pd 3 Mn, H and Pd 3 Fe is the Ll 2 structure, not the non- 
magnetically predicted DO22 structure. Thus, while the 
correlation between the calculated quantities N(Ep) ver- 
sus E(Ll2) — E(DC>22) holds, it leads to incorrect pre- 
dictions for the stability of Pd 3 Cr, Pd 3 Mn, and Pd 3 Fe. 
The generalized perturbation method calculations, ^] 
based on similar DOS arguments [diamond symbols in 
Fig. |) (a)] likewise predicts Pd 3 Cr, Pd 3 Mn, and Pd 3 Fe 
(and even Pd 3 Sc and Pd 3 Ti) to be stable in the DO22 
structures, in conflict with experiment. Q In this pa- 
per we explain this puzzle by noting that while a large 
N(Ep) indeed implies a destabilizing factor for the one- 
electron ( "band" ) energy, it also leads (in open-shell sys- 
tems) to spin-polarization and magnetic moment forma- 
tion which, in turn, is a stabilizing factor. Thus, despite 
their large N(Ep) in the LI2 structure (suggesting one- 
electron instability), Pt 3 Cr, Pd 3 Cr and Pd 3 Mn (nearly 
so for Pd 3 Fe) are correctly predicted to be more stable 
in this structure once spin-polarized total energy calcula- 
tions |[o| are done. Thus, magnetic ordering changes the 
predictions of NM total energy calculations and restores 
agreement with experiment. 

We have calculated the total energies of Pd 3 A, for X = 
Sc through Cu as well as Pt 3 Cr and Pt 3 Co in the LI2 
and DO22 structures using the LDA in both the spin- 
polarized and spin-unpolarized versions p"2[ of the full- 
potential LAPW method. p| In order to accurately 
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obtain the small energy difference between two fairly 
similar crystal structures, the calculations were carried 
out consistently using the same muffin-tin radii, Rmt, 
and basis set energy cutoffs, K max . The Brillouin Zone 
summations were done using the geometrically equiva- 
lent k-point sampling scheme" jl3|](a), in which 20 (40) 
k-points in the irreducible zone for the LI2 (EO22) struc- 
ture is mapped into the same 60 special k-points p3|(b) 
in the fee structure. We optimized the total energy as 
a function of volume, as well as the c/a ratio in the 
DO22 structure. The estimated LAPW error for the 
E(Ll2) — E{DC>22) energy difference is ~ 5 meV/atom, 
and the neglected zero-point energy difference between 
the two similar structures should be even smaller. At 
zero temperature, the absolute stability of a compound in 
a structure a with respect to phase separation is given 
by its formation enthalpy AH(o~) 
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FIG. 2. LDA calculated NM total DOS for Pd 3 X in the 
LI2 and DO22 structures. The inserts denote N(Ef). An 
asterisk denotes the more stable structure predicted by total 
energy calculations in the absence of magnetic ordering. Note 
that the more stable structure has a lower N(Ef). 
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AH a (a) = E a (a) - x A E A {V A ) - x B E B (V B ) (1) 

where Ea(Va) and Eb(Vb) are the energies of the con- 
stituents A and B in their respective ground states (e.g., 
for Pd 3 Ni, it is non-magnetic fee for Pd and ferromag- 
netic fee for Ni), E a (o~) is the energy of structure a, a — 
NM or FM denotes whether structure a is in the non- 
magnetic or ferromagnetic states. The relative stability 
of two different ordered structures is, 

SE a = E a (Ll 2 ) - E a (D0 22 ) , (2) 

while the magnetic stabilization energy of a given struc- 
ture a (LI2 or DO22) is 

SM(a) = E FM (a) - E NM (a) . (3) 

Table Q and Fig. ||(a) give SE^m and SEfm, while 
Fig. ||(b) shows the magnetic stabilization energies 
8M(Ll2) and SM(D022), and the local magnetic mo- 
ment fix on the X atom calculated numerically within 
the muffin-tin sphere (the value of nx is rather insensi- 
tive to the small change in muffin-tin radius). 



X = Sc Ti V Cr Mn Fe Co Ni Cu 
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FIG. 3. (a) LDA (present) and GPM (Ref. §) calcu- 
lated energy difference SE = E(L1 2 ) - E(DOi2) [Eq. (§)] 
for Pd 3 X compounds. Note the reversal of sign for SE due 
to spin-polarization for Pd 3 Cr, Pd 3 Mn, and Pd 3 Co (nearly 
so for Pd 3 Fe). Part (b) gives the magnetization energies 
5M(a) = E FM (a) - E FM (&) [Eq. @] for a = Ll 2 (empty 
circles) and a — DO22 (solid squares) structures and shows 
the local magnetic moment for the X atom (numbers above 
or below the symbols). 



Table || gives the formation enthalpies AH a (Ll2) and 
AH a (DC>22) for several systems. Figure || depicts the 
NM total DOS for Pd 3 X, where X = Sc->Cu, in the 
Ll 2 (left panel) and DO22 structures (right panel). One 
notices the following: 

(i) In contrast to the very similar DOS for the LI2 
and DO22 structures and the small energy difference 
E(L1 2 ) - E{D0 2 2) ~ 20 meV/atom predicted by the 
GPM, Q| one sees from Fig. |2] a marked difference of the 
DOS in the Ll 2 and DO22 structures, (a) While the DOS 
of the cubic LI2 structure resemble that of fee Pd, hav- 
ing three major peaks, the DOS of the DO22 structure 
are more smeared, reflecting a loss of cubic symmetry 
in the DO22 structure, (b) The DOS of the more stable 
Ll 2 structure for Pd 3 Sc and Pd 3 Ti shows a "pseudo-gap" 
near the Fermi level absent in the DO22 structure, (c) In 
the LI2 structure, the Fermi level of Pd 3 V, PdsCr, and 
Pd 3 Mn falls on a DOS peak (made mostly of d orbitals 
of the X atom), while in the DO22 structure the Fermi 
level falls on a relatively flat portion of the DOS. Indeed, 
these materials are more stable in the DO22 structure in 
a NM LDA description. As a result of these differences, 
the values of the N(Ep) (given in the inserts of Fig. ||) 
and its shape near the Fermi level are strikingly different 
for the LI2 and DO22 structures. 

(ii) The above noted trends in the NM N(E F ) in- 
duce a concomitant magnetic stabilization energy SM 
[Fig. |(b)] : A larger energy stabilization SM due to spin- 
polarization relates with a larger N(Ep) and with a larger 
localized magnetic moment \ix on the X atom [Fig. |^(b)]. 
For example, while the Pd 3 X compounds with X = Sc, 
Ti, and Cu are non-magnetic (so SM = 0), when X = 
Mn and Fe, one sees in Fig|(b) large energy lowering due 
to spin-polarization (SM ~ —200 meV/atom), and con- 
comitantly large magnetic moments of 3.5 (1b (X = Mn) 
and 3.1 [Lb [X = Fe) in both the LI2 and DO22 struc- 
tures (as a comparison, bec Fe has a magnetic moment 
value of only 2.2 (j,b)- Thus, spin-polarization induces 
a local magnetic moment on the "magnetic" 3d-atoms 
with large NM N(Ep) while lowering the total energy of 
the compound. In the case of Pt 3 Cr in the LI2 structure 
for which previous calculation exists, our calculated total 
magnetic moment in the unit cell (2.6 /j,b) agrees with a 

TABLE II. Non-magnetic (NM) and Ferromagnetic (FM) 
formation enthalpies, AH [in meV/atom, Eq. Aw], of some 
compounds. AH's are taken with respect to the NM fee Pd, 
Pt, FM fee Ni, FM fee Co, and anti-FM bec Cr, respectively. 
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previous calculation [Q of 2.6 /is and with the experi- 
mental value 1 15 1 of 2.5 [Ib- 

(iii) In the NM calculations, Pd 3 V and Pd 3 Cr have 
large DOS peaks at Ep in the Ll 2 structure, and are 
concomitantly less stable in this structure than in the low 
N(Ep) DO22 structure. However, as spin-polarization 
is introduced, the large DOS peak of the LI2 structure 
splits, so that Ep now resides in a low DOS region. This 
leads, simultaneously, to the formation of larger magnetic 
moments on V and Cr in the LI2 structure relative to the 
DO22 structure. This selective magnetization thus lowers 
the energy of the LI2 structure more significantly than 
in the DO22 structure [Fig. |(b)]. 

(iv) Table | shows that the more stable of the two struc- 
tures generally (and weakly) relates with a smaller value 
of the ferromagnetic DOS at the Fermi level Nfm(Ef). 
Thus, the trend of total energy stability with N(Ep) does 
exist, but for the spin-polarized quantities. 

(v) The magnetic stabilization energy SM, reverses the 
relative stability predicted by NM calculations in several 
cases: spin-polarization stabilizes the experimentally ob- 
served H Ll 2 structure of Pd 3 Cr, Pd 3 Mn, and Pt 3 Cr 
(nearly so for Pd 3 Fe) over the DO22 structure, while for 
Pd 3 Co, spin-polarization makes the DO22 structure more 
stable. The reversal of stability for Pd 3 Co can not be ob- 
served experimentally, since the calculated AH [Eq. ([!]) 
and Table |H| is positive, so Pd 3 Co (and similarly, Pd 3 Ni) 
is predicted to phase separate rather than to order, in 
accord with the observed phase-separation behaviors for 
Pd-Co and Pd-Ni. f§ 

(vi) Spin-polarization can stabilize ordering over 
phase-separation: jn| We find that in a NM descrip- 
tion Pt 3 Co has ATJnm > so it is predicted to phase 
separate, but that a strong spin-polarization effect sta- 
bilizes the ordered Ll 2 structure, leading to AiJpM < 
in agreement with the observed ordering behavior. |^] 
Similarly, spin-polarization stabilizes the experimentally 
observed [|17| LI2 structure of Pd 3 Cr (its magnetic be- 
havior has, however, not been experimentally examined). 
Hence, spin-polarization not only reverses the stability of 
the Ll 2 and DO22 structures for many compounds, but 
it also stabilizes an ordered (Ll 2 ) structure over phase- 
separation for Pd 3 Cr and Pt 3 Co. 

(vii) Interestingly, the AH for Pt 3 Cr are negative in 
both the NM and FM cases, but the spin-polarization 
effect gives a larger stability to the LI2 structure (ob- 
served experimentally) . [|| One also notices that the AH 
are lower in the Pt alloys than in the corresponding Pd 
alloys. The increased stability in Pt alloys has been ad- 
dressed previously in Ref. Jl6| . 

In summary, we have shown that a theoretically un- 
stable non-magnetic structures which involve magnetic 
atoms and possesses large N(Ep), may be stabilized by 
splitting the near Ep peak in the DOS and forming a 
local magnetic moment with a concomitant lowering of 
N(Ep) and the total energy. Therefore, theoretical stud- 



ies of the stability of compounds with a large value of a 
NM N(Ep) should be carefully tested for magnetic order- 
ing which can often change the predictions of the ground 
state crystal structure. 
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